Three Resonant Ultra-Cold Bosons: Off-Resonance Effects 
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We solve a finite range two-channel model for three resonant identical bosons. The model pro- 
vides a minimal description of the various magnetic Feshbach resonances in single species ultra-cold 
bosonic systems, including off-resonant scattering. We obtain important insights into the inter- 
pretation of seminal experiments: the three-body recombination rate measured in Sodium and the 
Efimov resonances observed in Caesium. This approach quantifies non universal effects appearing 
for a finite magnetic field detuning. 
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One of the main issues in current ultracold physics is to 
achieve highly correlated quantum gases. In these studies 
the magnetic Feshbach resonance is a crucial ingredient: 
it permits to tune the strength of the two-body inter- 
action measured by the s-wave scattering length to an 
arbitrarily large value, while the atomic density remains 
constant. However for dilute bosonic gases, such highly 
correlated states are very unstable as a consequence of 
three-body recombinations into deep molecular bound 
states [U, 13- Therefore, in the last decade an impres- 
sive effort has been made both theoretically and exper- 
imentally, for a deep understanding of three-body prop- 
erties in these systems. Universality concepts borrowed 
from Nuclear physics together with the specificity of ul- 
tracold atoms where an explicit energy scale separation 
occurs in scattering processes have led to unsuspected 
physical insights in this domain [3, For example, the 
so-called Efimov states have been observed for the first 
time in a resonant ultracold Caesium gas 01 • For scat- 
tering lengths much larger than the range of inter-atomic 
forces, Universal Theory Q permits a powerful analysis 
of three-body properties without any knowledge of the 
short range details of the real interatomic forces. For 
narrow resonances like the one observed in Ref.0], it is 
possible to fully determine the Efimov spectrum at res- 
onance using the effective range approach 1,0. De- 
spite their success, these approaches are not designed for 
a description of three-body properties at finite detun- 
ing, where off-resonant scattering effects come into play. 
Experiments on Efimov states in Ref. @ clearly exhibit 
some non universal behavior which cannot be taken into 
account by the Universal Theory. Moreover, the peak 
in three-body losses in Ref. [1] occurs for relatively large 
magnetic detuning where the effective range approach is 
inoperant. In this letter, we solve the three-boson prob- 
lem by using a two-channel model including the short 
range character of interatomic forces. The model de- 
scribes all the various type of magnetic Feshbach reso- 
nances: broad, narrow or in the neighborhood of a shape 
resonance. Results of the model compare quantitatively 
with experiments of Ref. 0] , highlighting the importance 



of off-resonant effects in narrow resonances. Concerning 
the experiments in Rcfs. 0,0, our results show that the 
observed violation of universality follows from the fact 
that the resonance in the atom-dimer scattering appears 
in a domain where the range of interatomic forces is not 
negligible as compared to the scattering length. 

We first introduce the finite range two-channel model. 
The two channels refer to the 'open' channel populated 
by atoms which are identical ultra-cold bosons of mass 
TO, and to the 'closed' channel where fundamentals enti- 
ties are couples of two tightly bound atoms that we here 
call 'molecules'. The model takes a simple form in the 
second quantized form, where the operator ak annihilates 
an atom in the open channel with momentum k, while 
annihilates a molecule of wavevcctor k in the closed chan- 
nel. Both at and bk obey standard bosonic commutation 
rules [ak,aj^,] = [^k,&k'] = (27r)'^(5(k — k'), corresponding 
to the choice (r|k) = exp(ik • r) for the plane wave state. 
Any other commutator vanishes. The Hamiltonian is 
similar to the ones introduced for two-component or fully 
polarized fermions in Refs. [lol [Tl|: 
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The first two terms in Eq. ([T]) are the kinetic op- 
erators in the open and closed channel respectively: 
Ek = h'^k'^ /{2m), and i?moi is the internal energy of the 
molecular state in the closed channel, defined with re- 
spect to the zero energy in the open channel. The mag- 
netic tunability of the interaction strength is due to the 
fact that -Binoi is an affine function of the magnetic field B 
with the slope where 5^ is the difference between the 
magnetic moments for an atomic pair in the open and 
closed channel. The second line in Eq. ([ij mimics the 
interatomic force in the open channel which is responsi- 
ble for the background scattering. In a real system it is 
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characterized by an attractive van der Waals tail, and the 
van der Waals constant Cg sets up the short range scale 
of this pairwise interaction: i?vdw — ^{mCQ/h'^y/'^ . In 
the present model, for simplicity we choose a separable 
interaction, with a Gaussian weight = exp(— fc^6^/2) 
imposing a cut-off at short distances: b = 0(i?vdw)- The 
last term in Eq. (jlj describes the coupling between the 
two channels and models the Feshbach resonance mech- 
anism (we choose A e M). The inter-channel coupling is 
also taken into account via the same Gaussian weight Xk 
as in the open channel interacting term. 

We now determine the different parameters of the 
model from measured two-body properties. For this pur- 
pose we solve the two-body scattering problem in the 
center of mass frame. For an incident plane wave of 
energy E = lT?k'^/m and momentum kg, the two-body 
state is a coherent superposition of one molecule plus 
two atoms: |*) = (/36jj -I- / Akaj^a^^)|0), where 



is the atomic wave function in the form: 



(27r)3j(k - ko) 
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and f{E) is the scattering amplitude. Remarkably, the 
s-wave scattering length (denoted by a) obtained in this 
two-channel model can be exactly identified with the fol- 
lowing expression which is known to be very accurate in 
the vicinity of a Feshbach resonance [l^ : 
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In Eq. ([3]), abg is the background scattering length, and 
AB is the width of the resonance located at the mag- 
netic field Bq. One finds a\^g = bgo^/^^ / {go — g^) , with 
gQ = —'in^/^h'^b/m, an energy detuning from resonance 
= £^moi — 2A^/ go + SfiAB where i/ = Sfi{B — Bo) and a 
resonance width AB = Snh^A^ahg/imgoSiJ,). This last 
relation implies that the energy width 6fiAB always has 
the same sign as the background scattering length abg. 
We checked that this property is indeed verified for the 
various resonances reported in Ref . [l2| • Using these pa- 
rameters, the scattering amplitude admits a simple ex- 
pression for E = —h?q^/m < (fco = iq and q > 0): 
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We use experimental or theoretical spectroscopic data 
on the two-body bound states as a way to choose a pre- 
cise value for b (in absence of such data, we arbitrar- 
ily set b = i?vdw)- In order to avoid any confusion the 
two-body bound states are denoted by 'dimers' and are 
distinct from the molecular state. Their binding energies 
E = —Edim = — ^^^dim/"^ < ^^'6 poles of f{E). In what 
follows, we briefly sketch their spectrum as a function of 
the energy detuning v. For < 0, there is a branch 
terminating at zero energy for v = 0^: this branch re- 
sults from the interchannel coupling and is denoted be- 
low as the Feshbach dimer's branch. Furthermore, in the 



case abg > by/ir, another branch exists for all possible 
values of the detuning. Away from the Feshbach reso- 
nance and for > 0, it results from the direct coupling 
in the open channel and we denote it as the 'background 
dimer' branch. For decreasing values of i> there is an 
avoided crossing between the two branches. We checked 
for several resonances that a choice of b of the order of 
-Rvdw permits to describe the lowest dimer's branch over 
a wide range of magnetic field detuning. In Fig. ([T]), we 
compare the results of the model with the experimental 
data of the resonance located at Bq ~ —11.7 G for Cae- 
sium {AB ~ 28.7 G, abg 1720ao, i?vdw ^ lOlao 
and Sfj, ~ 2.3|/xb| [l^). The system is in the vicinity 
of a shape resonance (abg ^ -Rvdw) and the spectrum 
displayed corresponds to the branch of the background 
dimer. We found very good agreement for b = 0.7i?vdW- 
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FIG. 1: Comparison between experimental dimer's energies 
in ^'^'^Cs obtained at Innsbruck (black losanges) Q and results 
of the two-channel model for b = 0.77ivdw (black line). Inset: 
detail of the deviation of the spectrum from the universal law 
Edim ~ /{ma^) (dashed line). 



We now turn to the central part of the present work 
by investigating the problem of three interacting bosons. 
The general three-body eigenstate of Eq. ([T]) of energy E 
in the total center of mass frame, is a coherent superpo- 
sition of three atoms plus one atom and one molecule: 
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In Eq.® ^K.k and /3k are respectively the atomic and 
the atom-molecule wavefunctions. For E > 0, the atomic 
wavefunction contains an eigenstate of the atomic kinetic 
operator ^^''k' ^^'^ G K, Ax.k can written as: 

^K.k = <|k + Axk/^K /(^^rci - 2ek + jO+), (6) 
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where E'rei = 'i^fcrei/™ = E — Sfi^ K"^ / (Am) is the relative 
energy between a pair of atoms and the third atom. For 
convenience, we introduced in Eq. (j6]) the effective atom- 
molecule wavefunction /3^: 

/3k = /?k X (S,ei -1^ + SfiAB)/{E,,,,i -v + S^iAB). (7) 

The three-boson problem for this model is solved by find- 
ing the solutions of the integral equation which is deduced 
from the stationary Schrodinger equation: 

4^fi2 /(Srcl) J (2^)3 Ek + CK + fK+k - E - iQ+ 

Remarkably all the two-body physics, contained in the 
scattering amplitude f{E), appears in the diagonal part 
of Eq. ([8|). Nearby a resonance, in the limit where 
\a\ — > oo, Eq. ^ converges asymptotically in the low- 
energy regime {i.e. for kb <^ I and Kb <^ 1) toward 
the so-called Skorniakov Ter-Martirosian equation of 
Ref. [l^ , while the high energy limit of the integral ker- 
nel acts as an ultra-violet cut-off and the Thomas collapse 
is avoided [l^. As the magnetic detuning |S — So| is in- 
creased, off-resonant effects come into play and the model 
can be used to quantify deviation from the universal the- 
ory 0, 01- To this end, we numerically solved Eq. ([8]) 
using two independent codes for several Feshbach reso- 
nances. We report in the following our results concerning 
the resonances studied experimentally in Refs. 
First, we investigate the spectrum of trimers made of 
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FIG. 2: Spectrum of bound states in "^Cs (So 11.7 G). 

Solid lines: trimers; dashed line: dimers; dotted horizontal 
line: "High energy" limit —h^/{mR^^y^); doted vertical lines: 
position of the Efimov resonances observed in the experiment 
(at 7.5 G and 25 G resp.); Crosses: thresholds computed with 
the two-channel model (at 10.3 G and 23.2 G resp.). 

Cffisium atoms for the Feshbach resonance at —11.7 G. 
Trimers are obtained by searching negative energy solu- 
tions (or E < —Edim if a shallow dimcr exists) of Eq. ([5]) 



where = 0. At resonance {v = 0), we recover the 
Efimov spectrum: En = Eoe~^'^"^'"' , (n e N) where Eq 
is the energy of the lowest Efimov state [of the order 
of-SV(TOi?2^w)] and So = 1.00624 In Fig. (0) we 

plot the spectrum as a function of the magnetic field. Our 
model predicts four important features: (i) The existence 
of two Efimov branches extending in the regions of posi- 
tive and negative B. Each branch is continuous through 
the formal limit \B\ = oo. Other Efimov branches are lo- 
cated in a small interval near Bq] (ii) The threshold of the 
first Efimov branch at 10.3 G observed in Ref.jH] at 7.8 G; 
(iii) The two trimer branches hit the background dimer 
branch in a non universal region where short range details 
of the interatomic forces are not negligible -see inset of 
Fig.([T])- and where our model gives qualitative informa- 
tions only; {iv) In experiments reported in Ref. the 
magnetic field was decreased from a large detuning and 
an atom-dimer resonance loss was found at ^ 25 G. The 
present model shows that the observed threshold corre- 
sponds to the second trimer branch. The disagreement 
between theoretical and experimental thresholds follows 
from item (iii). Moreover the asymptotic behavior of the 
wave function at distances r <C a is a crucial ingredient 
in Universal Theory and follows from taking the 'unitary 
approximation' of the scattering amplitude f ^ — 1/ifc at 
intermediate momentum 1/a <^ k <^ min(l/i?vdW: l/'"c), 
where is the effective range. However, at the ob- 
served threshold ^ 25 G, a/R^dw ^ 4 and from Eq. (|4|), 
o-ffc ^ 2. Hence, in this detuning region there is no clear 
separation of momentum scale and the 'unitary approxi- 
mation' for intermediate momentum is not correct. This 
fact gives an important insight in the deviation from pre- 
diction of Ref. [J] concerning the ratio between the scat- 
tering lengths at the two observed thresholds. 

In this last part, we study the recombination of three 
incoming atoms of vanishing total energy into one two- 
body bound state and one atom. This is the main process 
responsible for the short life time of a resonant BEG in 
a dipolar trap. Experimentally, it is measured via the 
atomic loss rate which is defined for N atoms trapped in 
a cubic box of size L by: 

7V--a,eciV(iV-l)(iV-2)/L6, (9) 

where arcc is the three-body recombination constant. 
The atomic loss rate plotted as a function of the mag- 
netic field detuning exhibits a large peak centered at the 
resonance. We limit our analysis to the standard regime 
where for a given detuning there exists at most one shal- 
low dimer only (obg is of the order of Rvdw or nega- 
tive). There arc two distinct regimes in the recombina- 
tion process: (a) In the regime of negative energy detun- 
ing (ly < 0) the formation rate of Feshbach dimers is the 
dominant loss mechanism. This dimer is an eigenstate of 
the model Hamiltonian in Eq. ([T]) and the recombination 
constant can be computed exactly; in what follows, it is 
denoted by af^^^^; (b) In the regime i/ > Q there is no shal- 
low dimer and the deep bound states populated in inelas- 
tic scattering processes are not described by our model 
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FIG. 3: Three-body recombination constant func- 
tion of the magnetic field B for tlie narrow Feshbach resonance 
in ■^■^Na at 907 G. The three different symbols correspond to 
different ramp speeds of the magnetic field across the reso- 
nance in Ref. Solid line: prediction of the two-channel 
model, dashed line: prediction of the effective range model in 
Ref. [y| . The vertical line at resonance is due to the existence 
of narrow oscillations of arcc also predicted in Refs. 0, Q]- 



Hamiltonian. In this case, we denote the recombination 
constant by afcc^- In both regimes the source term in 

Eq.® is A^^^ = (27r)6(5(K)(5(k) and corresponds to the 
incoming atomic wave function in the three-body scatter- 
ing process. In the first regime [case (a)], the outgoing 
dimcr and atom created by inelastic scattering have a 
relative momentum ^gdim- The dimer formation man- 
ifests itself as an outgoing wave in the atom-molecule 
wave function, and the wave function /3k has a pole at 
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iO'^; wc denote the residue of this pole by 



7. In order to evaluate ttrec*^, wc enclose the three incom- 
ing atoms in a fictitious box of arbitrary large size L and 
impose periodic boundary conditions. The wave function 
is then deduced from Eq. © with |^''^°'^) ~ \^)/L^/^. 
The rate of molecule formation is obtained in configura- 
tion space by computing the total flux of the probabil- 
ity current associated with the atom-molecule wave func- 
tion (relative particle of reduced mass 2m/ 3) through the 
box surface. Since the flying dimcr has a non vannishing 
probability (pciosed) to be in the closed channel, we di- 



vide this total flux by Pciosod thus obtaining the atomic 
loss rate N. By comparing this expression with the case 
= 3 in Eq. one obtains: 
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In the regime v > [case (&)], the recombination rate 
can be evaluated on a qualitative basis only. The idea 
is to compute in the fictitious box of size L the prob- 
ability V<: of finding the three atoms or the atom and 
the molecule in a volume of the order of i?vdw This 
calculation is performed in configuration space by using 
Eq. ^ . Since i?vdw gives the typical size of deep bound 
states, the loss rate is obtained by dimensional analysis 
with iV cx — __^J^>'^ , and finally the recombination con- 



stant is estimated by: aSijf'P = ^^nT'^ ■ We applied this 

formalism for the narrow resonance [i.e. SjiAB ^ ^^.^ ] 

in 23Na at So - 907 G Q (AS ~ 1 G, S^l 3.8\fiBl 
Obg — 63ao and R^dw — 44.5ao [111)- Fig- ® shows a 
dramatic agreement of our results (for b = Rvdw) with 
the experiments in Ref. 0]. Three-body properties in 
narrow resonances are usually described within the effec- 
tive range approximation by using the two parameters 
a and = h"^ / {mahgS fiAB) Q. As shown in Fig. (O, 
this latter approach gives reasonable results very close to 
the resonance only. In the zero range limit (6 — >■ 0), when 
the parameters a and R* are held fixed, and Obg = 0(6) 
vanishes (but at,g ^ by/n), the present model coincides 
with the effective range approach exactly. Therefore, in 
realistic situations where b and Cbg are both finite, this 
two-channel model quantifies consistently off-resonant ef- 
fects. We verified that arcc in Fig. ([3]) is not very sensible 
to the precise choice for b= 0(i?vdw)i showing that for 
a narrow resonance the present model gives quantitative 
results for the deviation from universality. 

To conclude, we presented a rather simple formalism 
to capture the main physical features in resonant three- 
bosons systems. The short range details of interatomic 
forces are described by one parameter only, so that we 
avoided the complexity of more detailed models [HI, [l^ . 
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